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^ ' Abstract 

o ■ 

Cosmological singularity and asymptotic behaviour of scale factor of generalized 
cosmological models are analyzed in respect of their structural stability. It is shown, 
CO ' that cosmological singularity is structurally unstable for the majority of models with 

barotropic perfect fluid with strong energy condition. Inclusion of A-term extends the 
set of structurally stable cosmological models. 
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1 Introduction 



Relativistic Friedmann-Robertson- Walker cosmological models (FRW-models) perform gen- 
erally accepted current theoretical basis for description of global structure and evolution of 
the Universe. Remarkable property of these models is their simplicity: in most situations 
' we deal with one scale factor and equation of state, which has to be specified by some addi- 

QQ , tional physical considerations. Observational data allow in principle to correct some basic 

■ parameters of the models, as well as to specify properties of matter at different stages of 

I evolution of the Universe and sometimes they imply reconsideration of Einstein gravity in 

the context of some generalized theories [Tl[21[3]. History of cosmology shows (and method- 
ological principle of eligibility prescribes) that correction of standard cosmological models is 
realized mainly by sequence of their, in some sense, "small" modifications. In this way we 
are necessarily faced with situation which is well known in qualitative differential equations 
theory [3]: some properties of original (non-modified) model may "survive" after "small" 
modification of the model, while others may disappear. In the former case the property is 
}J] ■ referred to as "rough" or structurally stable, in the latter case — as "thin" or structurally 

I unstable. As an example of such thin property of FRW-models we may consider cosmolog- 

ical singularity. By Hawking-Penrose theorem [5] which has been proved within Einstein 
GR, collapse is inevitably reached at some finite moment of time. However, there are many 
exact cosmological solutions, obtained in frame of generalized theories of gravitation (for 
example, involving scalar fields or non-riemannian objects [3 [7]), which describe evolution 
of the Universe without singularity or with singularities of qualitatively different kinds. By 
general considerations consecutive cosmological model should have structurally stable (with 
respect to possible modifications) basic properties. In the opposite case any "small" work 
modification of the theory will almost always lead to a change of cosmological paradigm. 

In present paper we perform analysis of structural stability of FRW-models. We restrict 
ourselves by investigations of structural perturbations induced by generalized /(i?)-theories 
of gravity and nonlinear equations of state. Both of these generalizations of standard cosmol- 
ogy have been actively investigated last decade in the context of new cosmological observa- 
tional data fS*, '9', IHl [HI [El [ISl [Hj . We investigate structural stability of the two important 
properties of FRW-models: cosmological singularity (S'-stability) and asymptotic behaviour 
of scale factor at t ^ cx) (A-stability). Latter property makes sense only for open and flat 
FRW-models. These restrictions allow to reduce the problem to investigation of asymptotic 
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solutions to second order differential equation with variable coefficients. In spite of quite 
different motivation of this work, our analysis in some aspects reproduces and supplements 
results obtained in [TTl [T31 [HI [111 [H] • In some important points we go to the conclusions 
made earlier in [161 [H] • 

In section[5]we outline class of standard models. In section[3]we describe those generalized 
models which will be used for investigation of structural stability of standard FRW-models. 
Basic equation on structural perturbations and its analysis are performed in section [31 
Section [5l is devoted to investigation of structural stability of standard FRW-models with A- 
term. Here we apparently illustrate the idea of "structural security" of cosmological model. 
In Conclusion we summarize and discuss all obtained results. 



2 Standard FRW-models 

Standard homogeneous isotropic cosmological models can be described by metric of the 
following kind: 

g = N'^{t)dt (g)dt- a^{t){dr (g) dr + {dd ® dO + sm"^ d(p (g) dif)), (1) 

where N{t) is factor, defining gauge of cosmological time, a{t) is scale factor, describing 
cosmological evolution of space lengths, k is curvature constant parameter, r is radial coor- 
dinate, 9 and are angle coordinates. The parameter k has dimension of inverse length and 
can be real, imaginary or null. In the first case we deal with closed cosmological models, 
in the second case — with open cosmological models, in the third case — with flat ones. 
In what follows all formulas will have universal sense for all possible values of parameter k. 
The case fc = sometimes will require rather simple limit procedure: k ~^ 0. 

Standard calculations give the following nonzero components of mixed Ricci tensor: 

□ 3 f Nd d\ 1 2 3 I ( Nd d P 
^'=N^[N^-a)' ^i^^2 = i?3-^|^^---2^j-2-^- (2) 

Scalar curvature is expressed by the formula^: 




■ ? 

a a 



? i ? 
a a'^ j 



7-^1-4- (3) 



It is commonly accepted, that source of gravity in cosmology is barotropic perfect fluid 
with energy-momentum tensor of the kind: 

T ^{e+p)u®u-pg, (4) 

where u = Ndt — 1-form of 4-velocity of matter, defining comoving reference frame, e and 
p — scalars of energy density and pressure of the matter, connected with each other by 
barotropic equation of state: 

p = as, where a = const. (5) 

Apart from perfect fluid for the more generality of our formulas we include into the standard 
models A— term (which, excepting section [5l we shall assume zero). Einstein equations: 



G$^R'^0--S^pR = KT^+AS'^0, (6) 



^Hereafter Greek indexes a,/3,7 denote space-time components of 4-tensors, Latin indexes i,j,k — their 
space components. 
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derived from Einstein-Hilbert variational principle for the action (S^n — action for matter, 
K = SttG/c^ — Einstein gravitational constant, * — standard dual conjugation): 

S[g] = -^l ^iR + 2A) + S^, (7) 

take the form of basic dynamical equations: 



3d2 .11 2Nd „a dM k 
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For the simplicity sake we accept conformal gauge of cosmological metric ([T]), which is 
defined by condition: N — a. Also we go to units, where k — 1 and c — 1. For chosen gauge 
formulas ©-([S]) take the form: 



and Einstein equations ([8]) are expressed by the formulas: 

With using ^TU\i it is easily to check validity of the relation: 

£^-3{p + s)-, (11) 
a 

having sense of energy conservation in adiabatically evolving Universe. For equation of state 
([5]) relation pT|) leads to the integral: 

where c — integration constant. Combining (fT2|) and first equation in (fTO|) . we go to the 
useful formula: 



Aa2 ,„ ^^^^ 
a^aM + — -P. (13) 

Differentiating this formula and expressing derivative d again through this formula, we go 
to another useful formula: 

„ ^ c(1^3a)^_3„_^2^^2^^3_ 
6 3 

By the same manner one can obtain the following useful formulas, which will be used in 
next sections: 



i?=(3a-l)ca-3(i+")-4A; j^o ^ ^(1 + 3a) ^-3(i+a) _ ^ „-3(i+a) _ 

(15) 



R = -3c(l + a)(3a - l)a-^^^+''\h^ + ^ - (16) 

V o o 

R = -3c(l + a)(3a - l)a-3(i+")(3(l + a)^:^ _ c(7 + 9q) ^„i_3„ _ 1 ^ 2)A). (17) 

6 3 
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Equation under A = can be integrated in elementary functionf0. The result has 
the form: 

r (^/H73sin[fc(l + 3a)t/2]/fc)2/(i+3"), a -1/3; 

«W = <^ (18) 

[ oq exp[±-\/c/3 — k'^t], a = — 1/3. 

Simple analysis of the integral (fT5|) shows, that in physical region (a > — 1) A— term plays 
no role in vicinity of cosmological singularity, but it dominates under a — > oo for flat and 
open cosmological models. 

At the end of this section we perform nonzero components of second covariant derivatives 
(acting on scalar functions) V = V^V^ in FRW-models: 

(fad d d _ I d — a/i 1 f d^ 2d d 

(19) 

3 Generalized cosmological models 

In present paper we consider generalized cosmological models of the following two kinds: 

1. Models with nonlinear equation of state; 

2. Models with nonlinear on curvature action for gravity. 

3.1 Nonlinear equation of state 

Lets consider barotropic fluid with equation of state of the following kind: 

p = ae{l + ae''), (20) 

where a and b are parameters, responsible for nonlinear properties of the matter. Under 
a —I- equation ((20)) goef[f| to the standard equation ([5]) . Physical nature of the nonlinearity 
can be concerned with the more complicated mechanisms of interaction of matter particles, 
quantum effects, unknown fundamental interactions, etc. 

If energy density of such nonlinear matter has structural perturbation 6e, then its pres- 
sure will have structural perturbation of the kind: 

Sp = aSe + aae''+^ . (21) 

In this expression nonlinear term is not varied, since we assume, that it already has first 
order of smallness (with respect to a). 



3.2 Nonlinear models of gravity 

Lets consider general scheme of nonlinear theories of gravity in cosmological context. For 
space-time of FRW-models there is the only scalar invariant — scalar curvature R. So, gen- 
eral kind of action for nonlinear cosmological models, generalizing Einstein-Hilbert action, 
can be expressed by the formula: 

S, = ~IJ (22) 

■^Though integral for world time: t = Jg a(t) dt can be expressed in elementary functions only for some 
particular values of parameters a and k. 

^For the case of dust-like matter a = it is necessary to use the following limit procedure in (I20I I: a 0, 
(T — ^ oo, aa < oo. 
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where f{R) — arbitrary function of scalar curvature. Standard variational procedure[f| leads 
to the following expression for variational derivative of (I22p with respect to metric: 



^ 4^"/3--i(K/3-VaV^+5a/3n]/'(i?)-La/3/W). (23) 



where Q is generalized Einstein tensor. Let f{R) can be expressed through the generalized 
Taylor decomposition: 

f{R) = 2A + R+ J2 KR', (24) 

where 2A = Aq, Ai = 1 and s — any real number, excepting and 1. Choice of the Ai 
corresponds to our main problem: to investigate structural stability of cosmological models 
in vicinity of standard FRW-models. Choice of the Aq will allow us to investigate structural 
stability of A-term cosmology in frame of general scheme. So, in space of nonlinear theories 
of gravity we are interested by those theories, which are characterized by the following 
infinitely dimensional vector of parameters: (. . . As_j, Aq, 1, Ag^, . . . ), when all As^ 
under i ^ 0, 1. In view of this aim it is expediently to decompose (with using (|24p ) tensor 
Q on standard Einstein tensor ([6]) and "perturbing" part, which can be interpreted (up to 
a constant) as energy-momentum tensor T" of nonlinear self-action: 

Gap = Gap + AGap = G„/3 + [kaP - VaV/3 + 9apn]F'{R) - ^ 5a/3 F{R), (25) 

where 

F(^) =2A+ J2 ^sR^ 

and circle over values relates them to unperturbed cosmological metric of those FRW-model, 
whose structural stability we are studying. 

Typical s-th member of this nonlinear perturbation in right-hand side of perturbed Ein- 
stein equation will have the following kind: 



-{AGsT^ = {TfTp = -s\s{Rt - V"V^ +5$ □) R^-' +^5'IW. (26) 



With using (fT2|) - (fT7)) and p9)) after some algebra, we go to the following expressions for 
non-zero components of nonlinear self-action tensor (there is no summation over repeating 
indexes!): 

{Tf % = A,(3a - l)*"ic*a-3.(i+a)(^^) ^ Sfai+3"), (27) 

where 

A0 = -(3a(2s2-3s+l) + 6s2-7s-l); B'i = s{s - + a); (28) 



Al^Al^Al^~~{s-^sa~-l)A'l- ^ ^ Bl = -{s + sa - a - ^)B'l (29) 
Note, that all formulas are valid for any real value of parameter s. 

^Under varying Lagrange density f{R)\/—g it is helpful to use Palatini's identity: SRajs = ^"[^^"a/S ~ 

^Note, that all As have different dimensions: [As] = . In spite of conditional sense of "smallness" for 

any dimensional value, our analysis in vicinity of the point (. . . , 0, 1, 0, . . . ) in parametric space is absolutely 
correct. 
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4 Equations for structural perturbations 



In order to derive linear differential equation for structural perturbations it is necessary 
to extract null order in the left-hand (geometrical) side of structurally perturbed Einstein 
equations: 

G{a{t))=T" + T"'^\ (30) 

where we consider now Einstein tensor as nonlinear differential operator G, acting on space 
of scale factors. Introducing structural perturbation of scale factor: a{t) — > a{t) + 5{t), we 
go to the formula: 

G(a(t) + S{t)) = G(a(i)) + SG{5{t), ait)) + 0(6), (31) 

where 

SG{S{t), a{t)) = ^G{a{t) + eS{t))U=o (32) 
de 

— " differential of Einstein operator" , which is linear differential operator on its first argu- 
ment. Direct calculations with using formulas PH)) . ([5^ and (fT^ - P^ lead to the following 
nonzero components of 6G: 

SG'^oi^,a) = (Vca-i"3"/3- fc2(5-f (/c^ _ 2ca^^-^'' /3)6^ ; (33) 



SGl {6, a)^^6~ A /3 - k^S + ""^^ ^ a~''-^"6. (34) 

Substituting it into equations (|30p and keeping in mind, that in null order on S these equa- 
tions lead to Einstein equations PH)) for unperturbed FRW-model which are satisfied iden- 
tically, we go to the following linearized system of equations on structural perturbations: 

a-^ a-^ 3 
Using formula (I2ip and excluding from (I35|) - ()36l) structural perturbation 5e, we obtain unique 
basic equation: 

25 + 2(3a - l)v/ca-i-3"/3 - k^6 + [Qak^ + c(l - 3a)a"i"3"/3),5 = (37) 

a3(a(T")0-f (r")}-afT£''+i). 
Going to the new variable x — a(t), equation (j37p can be transformed to the form: 

2x^icx-^-^"/3 - k^)6" + {6- x5'){6ak^ + c(l - 3a)x-^-^"/3) = (38) 

x3(a(r")° -I- (T")}) - a(7c''+ix"3("+i)(f+i)+3. 

Since sources of structural perturbations in first order act independently from each other, 
we can analyze structural stability of FRW-models by means of equation ([38| separately for 
every source. Limit correspondence considerations suggest that structural perturbations of 
scale factor must be determined by those particular solutions of differential equation ((38|) 
which vanish when source of perturbations vanished. 

^So, we omit general solutions to homogeneous equation (1381 (without its right-hand side). It is interesting 
to note, that this solutions define class of infinitesimal fantom conformal perturbations, which up to high 
orders of smallness transform apparently different metrics, originated from the same matter source, to each 
other. 
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4.1 Structural stability of cosmological models within nonlinear 
gravity 

Let consider structural perturbations of /(i?)-cosmological models. Perturbation from s-th 
member of Taylor decomposition in our notations takes the form: 

^'(«(m")o + iT^l) = A,(3a - l)^-ic^x-3«(i+")+3(A, + B,x'+^''), (39) 

where 

As = aA° + Al ^ {a + - s)^^; Bs = aB"^ + BI ^ {2a + ^ - s - sa)B°. 

Looking at equation (|38|) with right-hand side ([39| together with ((28|)-([29ll. one can imme- 
diately conclude: all standard cosmological models with isotropic radiation {a — 1/3) and 
vacuum-like matter (a = —1) are structurally stable in any order on curvature, since this 
concrete values of a imply vanishing of sources of structural perturbations. 
In general case let introduce the following notations: 

. Six) . p _ 3fc2 



3As(3a - l)«-ic^-i(l -I- a){l - s)x ' 



c 



General condition of structural stability is boundedness of relative perturbation ^ (x) for any 
allowed x. Note, that in this section we accept, that a ^ —1 and a ^ 1/3 in view of above 
mentioned specific properties of FRW-models with such parameters. Also we put anywher^ 

In order to analyze iS-stability let consider asymptotic kind of equation (j38p together with 
its right-hand side ([55]) under x — > and under a > —1/3, when cosmological singularity 
does take place (solution p8)) ). For general values of parameters this asymptotic kind is 
expressed by the following equation: 

2x^ix^ix))" + {l-3a)x{^{x)~{x^{x)y) = i(3a(2s2-3s + l) + 6s2-7s-l)a;-3^(i+")+4+3". 

(40) 

Its particular solution, which we are interested, has the form: 

^(x)U.o = .^fl-^^^^^^^^f "^p^^. x--^-^^^^-^-^ (41) 
6(6s -|- Qsa — 7 — 9a)(l -I- a)(s — 1) 

So, condition of S'-stability is positiveness of index: — 3(s— 1)(1 + q;) > 0, that under l + 3a > 
takes place, when s < 1. In other words, almost all nonlinear analytic generalizations of 
standard cosmology possess structural instability of singularity. Note, that non-analytic 
nonlinear models (for example, containing in lagrangians negative powers of R) may be 
iS-stable. 

The words "almost all" are necessary to separate A-term cosmology (s = 0) and some 
particular cases, which are iS-stable. S'-stability of the particular case, when coefficient in 
([1T|) vanishes: 

6as'^ - 9sa + 3a + Gs^ - 7s - 1 = 0, (42) 
corresponds to the asymptotic equation (I40|) with other right-hand side: 



2x^{x^{x))" + (1 - 3a)x{^{x) - (x^ix))') = 3Ps(2a + 1 - s(l + a))a;-3^(i+")+^+<'". (43) 



^Situation with s = 1 has no particular interest, since, in fact, it can be reduced to redefinition of Einstein 
gravitational constant and so it can not infiuence on stability. 
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Using relation between parameters: 



a = -, — r , (44) 

3(2s2 - 3s + 1) ' ^ ' 

which foUows from ([32]), we obtain the following particular solution to equation (|43p : 

Omitting the situations with s < 1, when almost all models are S'-stable, we obtain two new 
classes of S'-stable models. The first one is obtained under s E (l;s*], that corresponds to 
monotonic variation of a G (oo; a,]. Values of boundaries of the intervals are: 

2 3(19-5^/T3) ^ ' 



Second class is obtained for flat models fc = 0. The interval of nontrivial 5-stability (i.e., 
interval, where a > —1/3) is described by the formula ([44]) for s G (1; (1 + y/3)/2). 

Finally, the last particular case corresponds to the following relation between parameters: 

6s — 7 7 + 9a 

(6s + 6sa-7-9a) = 0^a = -— ^ s = — — ^, (47) 

3(2s — 3) 6(1 + 0;) 

when denominator in (I4ip vanishes. Asymptotic particular solution in this case has the 
form: 

^(x)|_o = (l-3«)((l + 3«)M") + 2)^_a+3.)/2^ (48) 

^1x^0 2(l + 3a)2 ^ ' 

Under a > —1/3 (1 < s < 3/2) power is negative and the models are S-unstable, while 
under s < 1 there is no singularity and the models are S'-stable. 

So far the case a = —1/3 has remained out of our consideration. Particular solution to 
the equation ([55]) with right-hand side ([55)) . taken under a — —1/3, has the form: 

-2s2 + 2s + l-2fc2s + 2fc2s2 

- — sik^-Dis^ir — " ' ^''^ 

that leads to the conclusion: under s > 1 generalized nonlinear models with equation of state: 
p — — e/3 are S -unstable and A-stable and vise verse: under s < 1 such models are S -stable 
and A-unstable. Particular situation arises when coefficient in expression (j49p vanishes. In 
this case the model is structurally stable in all senses, at the least, in linear approximation. 
In closed models this situation is described by relation: 



fc2 - 1 ± Vfc4 _ 4fc2 + 3 

— w^) — ■ ^ ^ 

Since solution (fT5)) must be real, parameter k lies in interval (0; 1). Parameter s (on both 
branches of square root) belongs to the union (-oo; (1 - V3)/2] U [(1 + V3)/2; c5o). For open 
models relation between s and k can be obtained from (|50|) by the substitution k^ —k^. 
When k varies from to oo parameter s monotonically falls from (1 -I- •\/3)/2 to 1 on upper 
branch of square root and monotonically increases from (1 — \/3)/2 to on lower one. The 
value = 1 is singular for equation ()38|) . since under a ~ —1/3 and k^ — 1 this equation 
become algebraic: a;~2«+3 = Q. It corresponds to the limit s — > oo, which is reached on the 
upper branch of solution ([50)1 . In this degenerate case cosmological models can be both 
structurally stable, and structurally unstable, since basic equation (|38p is satisfied by any 
function 6{x). All results are summarized in the diagram of S'-stability (figH]). 
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Figure 1: Diagram of S'-stability of FRW-models. Inclined dash points region of 5- instability of 
FRW-models. Boundary of this region is double line. At this boundary FRW-models can be both 
S'-stable (vertical line a = 1/3, horizontal line s — 1 and part of rational dependency (|44p . beginning 
at the point {st;a,} (formula (I46|) ) and conditionally S'-stable (vertical line (with horizontal dash) 
a = —1/3 with condition (|50p : upper part of this line correspondes to closed models, lower — to the 
open ones; and part (with vertical dash) of dependency (|44p . lying between s* and s,, = (1 -|- y/3)/2 
with condition ~ 0)- 



In order to analyze A-stability of FRW-models it is necessary to go to asymptotic kind of 
equation ([38]) and of its right-hand side ((39|) under a; ^ oo. Note, that energetic condition: 
1 + 3a > is unnecessary now. 

In first, lets note, that under 1 -I- 3 a < asymptotic kind of equation for structural 
perturbations exactly coincides with ([^0]) . while its particular solution coincides with ([1T|) . 
which one should consider now as solution under x oo. Analysis of power sign leads to 
the conclusion, that under —1 < a < —1/3 almost all FRW-models are A-stable when s > 1 
and are A-instable when s < 1, and under a < —1 vise verse: almost all FRW-models are 
A-instable when s > 1 and A-stable when s < 1. Lets go to the particular situations. The 
case, when coefficient in (PT|) vanishes is described by the formulas ([1^ - (|^5|) . which should 
be considered as asymptotic expressions under x — s- oo and a < —1/3. Additional A-stable 
branch of FRW-models arises here when s £ (— oo; s*] U (2; oo), where s* = — (1 + VTS) /2 w 
-2.30. Here on the first interval a € (-l;a*], where a* = -(4 -I- VT3)/9 w -0.85 and a 
monotonically increases. On the second interval a decreases from —1 to its minimal value 
Q!min = (4V6 — 13)/3 ~ —1.07, which is reached under Smin = 2 + V6/2 w 3.22, and then 
monotonically increases to the value —1. There is no new A-stable models in case of vanishing 
of denominator in ([1T|) . 

Now lets go to the case 1 -I- 3a > 0. Asymptotic kind of equation for structural pertur- 
bations under takes the form: 

-2x^{x^{x))" + 6axi^{x) - (x^ix))') = 3s(2a + 1 - s(l + a))a;-3^(i+")+4+3a_ (5^) 

o 

Its particular solution: 

-3(s-1)(1+q) 

6(s - 1)(1 + a) 

So, under a > —1/3 almost all non-flat FRW-models are A-stable under s > 1 and are 
A-unstable under s < 1. In particular case s — Q equation (|5ip will have another right-hand 
side and another particular solution: 

(:( \\ — 3a — 1 2 
24(7^(1 -I- a) 
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that implies A-instability of all such models. 

Finally, for the case of flat models equation for structural perturbation become identical 
to (PD|) and its particular solution become identical to (|¥T|) . So, the models are A-stable 
under s > 1. The relation ((44)) defines additional branch of A-stable models in the region 
a > -1/3 under s** < s < 1, where s** = (1 - \/3)/2 « -0.37, while relation (gT]) defines 
additional branch of A- unstable models under 3/2 > s > 1. 

All results are summarized in diagram of A-stability (fig. [5]). 
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Figure 2: Diagram of A-stability of FRW- models. Inclined dash points region of A- instability 
of FRW-models. The boundary of this region is double line. At this boundary FRW-models can 
be ^-stable (vertical line a = 1/3, horizontal line s = \ and parts of rational dependency (|44p 
with end at the point {a*; s*} and the beginnings at the points { — 1/3; s**} and { — 1; 2} (formula 
(|46[)). conditionally A-stable (vertical line (with horizontal dash) a = —1/3 under (|50[) : whole 
performed part of this line corresponds to open models) or A-unstable (part (with inclined dash) of 
the dependency (|47p . beginning at the point {—1/3; 1}; empty point on this branch has coordinates 
{1/3; 5/4} and corresponds to intersection of (|47|) with vertical line q = 1/3, where all models are 
stable), yl-stable branch, which asymptotically tends to line s = 1/2, corresponds to = 0. 



4.2 Structural stability of FRW-models with nonlinear matter 

Let consider now the equation ([55]) with right-hand side of kind 2;-3(Q+i)(b+i)+3^ responsible 
for perturbation by nonlinear properties of matter. In the limits a; — s- under 1 + 3q; > 
this equation can be reduced to the form: 

2x^{xi{x))" + (1 - 'ia)x{i{x) - {xi{x))') = a;-3(i+«)(b+i)+4+3a^ (53^ 
where (,{x) — —iix^joLG&^^x. Its particular solution: 

^-3h(l+Q) 

^^"^^'^^^ " 36(-3a + 6a6-l + 66)(l + a)- ^^""^ 

Simple analysis of index shows, that almost all FRW-models with nonlinear matter are S- 
stable under 6 < 0, and are S -unstable under b > (see also interesting analysis in pi 121])). 
Particular solution to the asymptotic equation (|53p under 6 = has the form: 
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We see, that such models are S'-unstable. Particular solution to the equation ([55)1 under 
6 = (1 + 3a)/6(l + a), when denominator of (j54p vanishes, takes the form: 



(3a^l)hix + 2 (i+3„)/2 



(56) 



and we see, that such models are 5-instable too. Results of our analysis are shown in 
diagram (fig. [S]). 
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Figure 3: Diagram of S-stability of FRW-models with generalized nonlinear equation of state. 
Inclined dash shows region of S-instability of FRW-models. Boundary belongs to the region. 

In case (l + Sa) < analysis of A-stability is reduced to previous formulas of this section, 
if one will consider their behaviour under x — > oo. So, from the expression ()54p it follows, 
that almost all FRW-models with —1 < a < —1/3 are A-stable under b > and are A- 
unstable under b < 0; almost all FRW-models with a < —1 are A-stable under b < and 
are A-unstable under b > 0. Formulas allow us to conclude, that in all particular 

cases, which have been considered above, the models are A-unstable. 

For the case 1 -I- 3a > we obtain the following asymptotic form of equation 



-2Px^x^{x))" + 6aPx{C{x) - (x^ix))') = a;-3(i+a)(6+i)+3^ (57) 
Its particular solution has the form: 

rj^—l—3ah—3a — 3b 

? (x) U^oo = - g^^^ ^ a)~k^3ab + 3a + 35 + 1) ■ ^^^^ 

Analysis of index shows, that almost all FRW-models with strong energodominancy condition 
are A-stable under 

6< -(l + 3a)/3(l + a) (59) 

and are A-unstable in case of opposite inequality. 

Further analysis of particular cases of the models with strong energodominancy condition 
reveals that: 

• under 6 = models are A-stable; 

• under b = —(1 + 3a)/3(l + a) models are A-unstable; 

• under k = Q models are A-stable under 5 > 0, while under 6 < models are A-unstable; 

• under k = and a = —(66 — l)/3(26 — 1) models arc A-stable; 

• under a = —1/3 models are A-stable and S'-unstable when 6 > 0, A-unstable and 
iS-stable when 6 < 0. Under 6 = models are unstable in all senses. 

• under a = — 1 models are A-unstable under any value of 6. 
All results are summarized in diagram of A-stability (fig. |4]). 
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Figure 4: Diagram of j4-stability of FRW-models witli generalized nonlinear equation of state. 
Inclined dash shows region, where the models are yl-unstable. Boundary of this region belongs to 
the region. Main part of this region on lower half-plane is bounded by the curve (|59|l . Half- line 
6 = 0, Q > — 1/3 corresponds to flat >1- unstable FRW-models. 

5 Stability of standard models with A-term 

So far we have considered non-perturbed FRW-models with A = 0. One may hope, that 
extension of class of non-perturbed models will lead to the extension of their stability region. 
Let us illustrate this idea by example of non-perturbed FRW-models with non-zero A-term. 
For our purposes it will be sufficient to restrict ourselves by FRW-models with dust-like 
matter (a = 0) for generalized nonlinear theories of gravity and by the FRW-models with 
isotropic radiation (a = 1/3) for generalized models with non- linear matter. By note, 
made after formula (HH]), inclusion of A-term does not change behaviour of FRW-models 
near singularity, that may be checked directly by our method. The diagram [5] implies, 
that generalized FRW-models with high orders of scalar curvature under a = (and 
under A = 0) are A-stable under s £ [l;oo) \ {7/6} for all values of A:^ < and under 
s = (7 — \/73)/12 « —0.13 for k = 0. The diagram |4] implies, that generalized FRW-models 
with isotropic radiation, perturbed by non-linearity of state equation (under A = 0), are 
A-stable under be [— l/2;oo) for all values oik^ < 0, excepting the value 6 = under fc = 0. 

Standard calculations with using formulas (fTO|) - pTl) . taken under A 7^ 0, lead to the 
following kind of equation on structural perturbations for models with non-linear gravity 
under a = in the limit x 00 : 

where we have omitted constant factor in the right-hand side. Its particular solution: ^ = 
i5/a; ~ const, that means unconditional A-stability of all FRW-models with A-term and dust- 
like matter. 

For FRW-models with A-term and nonlinear equation of state we obtain equation for 
structural perturbation of the following kind: 

x^S" + 2xH' - Qx'^5 ~ x-^-^\ 

Its particular solution: ^(x) ~ x~^^^^^\ So, we conclude, that due to inclusion of A-term, 
region of A-stability of FRW-models with isotropic radiation become larger: it includes now 
additional segment of b € [—1; —1/2] and the value b — 0. This conclusion is closely related 
to results of papers [l6j[T7], where generic structural stability of ACDM cosmological models 
has been established by analysis of topological properties of phase portraits of cosmological 
models. 

Complete analysis of structural stability of FRW-models with A-tcrm wc defer for future 
pubhcations. 
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6 Conclusion 



We have shown, that the question of structural stabihty of standard cosmological FRW- 
models is of non-trivial nature and can be investigated in general manner in frame of out- 
lined classes of generalized cosmological models, including non-linear gravity and non-linear 
properties of matter. We have restricted ourselves to investigations of structural stability 
near the cosmological singularity under a{t) —> (5'-stability) and under a{t) oo (A- 
stability). Note, that in terms of relative structural perturbation 5-stability condition takes 
the form: 

lim - < oo. (60) 

a^O a 

So, S'-instability condition can be considered in the following two qualitatively different 
situations. As it can be seen from the content of sections |4] and [5l in almost all cases 
asymptotic behaviour of structural perturbation is described by power function of a: (5 ~ a^. 
The condition ([60|) leads then to the inequality p > I, providing S'-stability. It is obvious, 
that under < P < 1 singularity is conserved, but behaviour of scale factor near singularity 
changes. In case p < 0, singularity is destroyed. In our diagrams of S'-stability these 
situations are indistinguishable. In case of ^-instability we, obviously, always deal with 
change of asymptotic behaviour of scale factor. 

Union of all stability diagrams, which corresponds to simultaneous action of structural 
perturbations of all considered types, leads to the following conclusion: only the models with 
— I < a < —1/3, 6 > 0, s > 1 and the models with a < —1, 6 < 0, s < 1 together with some 
particular families of cosmological models, having zero measures on diagrams [TKH possess 
complete (in frame of chosen class of perturbations) structural stability. It, in particular, 
means that almost all FRW-models with strong energy condition are structurally unstable. 

In our opinion the question of structural stability in cosmology is of principal signifi- 
cance. Consequent account of new observational data and evolution of our theoretical tools 
of descriptions of nature lead to unavoidable theoretical modifications of basic cosmologi- 
cal equations. If we don't speak about scientific revolution and about change of paradigm, 
such modifications will be looked at as, in some sense, small corrections in lagrangians or in 
equations. Our analysis shows that in spite of " smallness" , this corrections may lead to the 
total deleting of one of the properties of original theory, to significant modifications of oth- 
ers properties and to insignificant quantitative variations of the third ones. In the first two 
cases the properties of original models are "thin" with respect to the modifications of the 
theory and so, in the context of generalized models, they are "accidental" . Our investigation 
reveals, that cosmological singularity is such "accidental" property of Einstein cosmology. 
Any non- linear corrections of kind \W under s > 2 in gravitational lagrangian will destroy 
cosmological singularity, whatever small constant A is chosen. In fact, cosmological singu- 
larity is not experimentally observed. It is desirable to provide stability of those properties 
of cosmological models which have robust experimental (observational) support. In spite of 
infinity (and even uncountability) of a set of possible generalizations of standard cosmology, 
the set of its viable and reasonable (from the viewpoint of present time) generalizations is 
quite foreseeable and can be analyzed from general positions [1] . 

Our example with A-term shows that among the set of allowed models one can try to 
find relatively simple one, which possesses sufficient "structural security" and conserve some 
important properties under its possible further modifications. Note, that considerations 
of structural stability can be viewed as purely theoretical argument for inclusion of A- 
term into Einstein theory of gravity. Lets note also that diagrams [2] and H throw light on 
the following question: what class of generalized cosmological models can be relevant to 
accelerated in latter times universe? The answer is the following: for FRW models with 
strong energodominancy condition (1 -|- 3a) > such models must be A-unstable! In other 
words relevant generalized models must have in their lagrangians terms of the kind i?* with 
s < 1 and (or) equation of state with nonlinearity index b < —(1 + 3q;)/3(1 + a). 
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Using technic of conform transformations [22 our considerations can be directly applied 
to the cosmological models within Einstein gravity with self-interacting scalar fields and 
within non-riemannian geometrical theories. This topic is far from the scope of present 
paper. 

In fact, the question of structural stability is not limited only to cosmology. It is valid 
in all situations when physics is local and can be described by sets of differential equations. 
By general considerations any local physical model can be "slightly" modified, such that 
any forgiven property of the model will be destroyed. Probably, in this situation some 
fundamental principles will be helpful to restrict the set of possible modifications of the 
model and to "stabilize" some of its important physical properties. 
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